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UPPER HANUE HETEROGENEITY: COMPARISONS OF REGIONS SOUTH OF AUSTRALIA 


WITH PHILIPPINE BASIN 


The goal of this study has been to attempt to identify the nature of 
h £ mass anomalies that ray occur beneath the regions of negative residual depth 
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positive gravitational attraction due to the excess density of the cold 
mantle woul d be expected, and might provide the compensation for the mass 
deficiency owing to the excess water depth. Because the topographic mass 
anomaly is closer than the upper ran tie mass anomaly, observed free -air 
anomalies would be slightly negative even for perfect isostasy. 



model studies of convection In the upper mantle also predict a positive 
correlation between residual elevation and free-air gravity anomalies. No 
obvious correlation with long wavelength free-air gravity anomalies has been 
observed south of Australia, and a region of slightly negative anomalies 
occurs in the West Philippine Basin. 

Under other funding, we have recently identified that for each spherical 
harmonic degree, the ratio of the gravity contribution to the geoid contri- 
bution Is a greater constraint on dimensions and depth of possible source 
masses than Is the wavelength associated with the zonal harmonics. This 
realization has led to a plausible decomposition of the earth's gravity 
field (Bowin, 1980, in prep.). This decomposition suggests that the mass 
anomalies due to topography at the core-n»ntle boundary are predominantly 
contained in the coefficients for the second and third harmonics, and that 
the harnonics 4 through 10 contain information on nnss anonul ies in the 
mantle -- largely below 300km depth. The combined contributions from degree 
4 through 10 yield a narrow belt of positive geoid anomalies coincident with 
most of the world's sites of plate convergence. That positive anomaly may 
result from nuss excesses associated with subducted cold lithosphere. By 
subtracting a degree 10 field from radar altimeter geoid observations, the 
contributions from mass anonul ies at the core-nan tie boundary, and from the 
deeper parts of subduction zones and of the mantle are largely removed. 

The residual geoid anonul ies are thus restricted to those resulting from 
mass anomalies in the outer 600km of the earth as would all mass anonul ies 
contributing to degrees 11 and higher on the basis of their gravity to 
geoid ratio. 


Vo better evaluate the specific regional spherical harironlc field to re- 
move In localized studies we plotted cummulative contribution curves for the 
Philippine Basin and the region south of Australia (Figure 1 and 2). From 
the cummulative contribution curves, a spherical harmonic degree 10 field 
was selected for the Philippine Basin for subtraction of a regional field 
best representing mass anomalies of the earth and the very broad wavelength 
anomaly presumed to be due to the subducted lithosphere slab at the Marianna 
arc/trench system. A spherical harmonic degree 13 field was selected for 
the region south of Australia. These spherical haruonic representations of 
a regional field were subtracted from GE0S-3 geoid data. Contoured maps of 
the resulting residual geoid maps are shown In Figures 3 and 4. 

The residual geoid anomaly maps show a positive anomaly associated with 
the residual depths south of Australia, and a negative anomaly associated 
with the West Philippine Basin. The residual geoid anomalies are consistent 
with the anomalous surface wave velocities, and with a thermal origin, but 
obviously not then with residual depth. The pattern of positive residual 
geoid anomaly In the region south of Australia, however, differs from the 
northeast trending negative residual depth anomaly as presented by Cochran 
and Talwani (1977). The site of the greatest residual depth anomaly 
(-1 ,000m) lies beneath the east-west trending residual geoid low Immediately 
south of the Australian shelf. Further examination of the data Is required. 

The residual geoid low in the Philippine Basin Is consistent with the 
area of negative free-air gravity anomalies seen In the map of Watts et al . 
(1978) and the area of deep water (depths greater than 5900m) In the West 
Philippine Basin (Mammerickx, 1976). Comparing the ratio of the regional 
free-air gravity anomaly to the residual geoid anomaly with that due to 


thin horizontal disks suggested a depth to the anomalous irass source of about 
200knt; that Is, considerably below the base of the lithosphere. However, 
nore rigorous model Ing, particularly one incorporating an upper mantle mass 
excess as isostatic compensation for excess water depths was needed. We 
pursued such an Investigation using theoretical admittance functions, and a 
preprint of the resulting paper is appended to this report. A summary of 
these results will be given at the end of this report. 

Under this grant we have recently succeeded in developing an interpola- 
tion scheme utilizing correlation functions that together with our geographg- 
1c organized data files allow estimation of geoid and gravity values at any 
selected location. The first results obtained using this technique are shown 
in Figures 5, 6 and 7: profiles of GEOS-3 geoid, free-air gravity anomalies 
from surface ship data, and bathymetry along Latitude 20°N were construc- 
ted. Computed GDI-9 geoid and gravity fields to degree 10 were then sub- 
tracted from the Interpolated values to obtain the profiles of residual geoid 
and gravity along 20°N latitude. We have subsequently learned how to 
better tune the Interpolation parameters and now can obtain better estimates 
(See Fig. 8 for example of improvement). This interpolation technique also 
provides error estimates for each value (Fig. 8). This capability Is very 
important for quantitative analysis of the geoid and gravity data. 

The ability explained above to retrieve good estimates of geoid and 
gravity anomalies at preselected locations and thus construct any desired 
profile or grid will be used to experiment in the contained Inversion of 
residuai geoid and gravity data to obtain estimates of the dimensions and 
depth of the mass anomalies causing the geoid and gravity anomalies. 


In summary, residual geoid anonal ies in association with negative resi- 
dual depth anomalies have been identified in the Philippine Basin (negative) 
and in the region south of Australia (positive and negative). In the latter 
region the geoid anomalies have eastward trends whereas the depth anomaly 
trends northeast. The surface wave velocity data should now be re-examined 
to see hay they may fit to this nay finding. The residual geoid anomaly is 
positive in the vicinity of the crest of the southeast Indian ridge where 
the anomalously fast surface wave v,.T *ci t> data seems best defined, and 
anomalously cool upper mantle could explain both. However, if this cool 
mantle were the cause of the excess water depths, negative geoid and gravity 
anonal ies would be expected. In the Philippine Basin, the negative residual 
geoid anomaly is consistent with a thermal cause for the slow surface wave 
velocity structure, but if so, then a warm (slow) upper mantle would not be 
an explanation for the excess water depth (negative residual depth) 
anomalies in the Philippine Basin. 

Theoretical gravity -topography and geo id -topography admittance functions 
indicate that for the Philippine Sea region high density mantle at about 60 
km depth can account for the nagnitudes of the gravity and residual geoid 
anomaly and the 1 km residual water depth anomaly. In that case, however, 
dense upper mantle material would be contrary to the anomalously slow 
surface wave velocity structure in the Philippine Sea as determined by 
Seek ins and Teng (1977), although it could be compatible with a normal 
surface wave velocity structure in the west Philippine Basin as determined 
by Sacks and Shione (1981). Alternatively, the negative residual depth 
anomaly may be compensated for by excess density in the uppermost mantle, 
but the residual geoid and regional free-air gravity anomalies and a slow 



surface wave velocity structure might result from low-density warm upper 
mantle material lying beneath the zone of high-density uppermost mantle* 
From a horizontal disk approximation, the depth of the low-density warm 
mantle is estlmted to be on the order of 200 km. 

Both the negative depth anomaly and the compensating mass excess In the 
uppermost mantle nay have developed in the Eocene as the lithosphere of the 
west Philippine basin formad. Heating of the deeper upper mantle which 
could cause slay surface wave velocities and the negative gravity and geoid 
anomalies may be a nuch younger phenomenon and presently be in progress. 
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Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 
Figure 7 


Figure 8 


Cumulative geoid anomaly curves for 15 sites In the Philippine 
Basin region. Computed from GEM-9 spherical harmonic 
coefficients. Numbers below lower right curve indicate 
harmonic degree. 

Cumulative geoid anomaly curves for 12 sites in the region 
south of Australia. Computed from GEM-9 spherical harmonic 
coefficients. Numbers below lower right curve indicate 
harmonic degree. 

Residual geoid anomalies of the Philippine Basin region. 
Obtained by removing a GEM-9 spherical harmonic degree 10 
field from GEOS-3 observations. Contour interval is 1 meter. 

Residual geoid anomalies of region south of Australia. 

Obtained by removing a GEM-9 spherical harmonic degree 13 
field from GEOS-3 observations. Contour interval is 1 meter. 

Profiles of geoid anomaly, free-air gravity anomaly, and 
bathymetry across the Philippine Basin along Latitude 20°N. 
Data obtained by interpolation using objective analysis 
techniques. 

Profiles regional geoid and gravity profiles computed from 
GEM-9 spherical harmonic coefficients through degree 10. Same 
location as for Figure 5. 

Profiles of residual geoid anomaly, residual free-air gravity 
anomaly, and bathymetry, obtained by subtracting regional 
fields (Fig. 6) from observed data (Fig. 5). 

Profile of geoid anomaly across the Philippine Basin along 
Latitude 20®N. Interpolated values as in Fig. 5 but with 1% 
noise assumed (top). (Middle) profile Is number of GEOS-3 
observations used in each interpolation. (Bottom) profile is 
predicted error In dimensionless values reflection 
distribution of data observations for each interpolation. 
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